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This paper extends an earlier analytical scattering matrix treatment of conductance and local- 
ization in coupled two- and three Anderson chain systems for weak disorder when evanescent states 
are present at the Fermi level. Such states exist typically when the interchain coupling exceeds 
the width of propagating energy bands associated with the various transverse eigenvalues of the 
coupled tight-binding systems. We calculate reflection- and transmission coefficients in cases where, 
besides propagating states, one or two evanescent states are available at the Fermi level for elastic 
scattering of electrons by the disordered systems. We observe important qualitative changes in these 
coefficients and in the related localization lengths due to ineffectiveness of the evanescent modes for 
transmission and reflection in the various scattering channels. In particular, the localization lengths 
are generally significantly larger than the values obtained when evanescent modes are absent. Effects 
associated with disorder mediated coupling between propagating and evanescent modes are shown 
to be suppressed by quantum interference effects, in lowest order for weak disorder. 

PACS numbers: 72.15.Rn, 73.21.Hb, 73.63.Nm 

I. INTRODUCTION 

In a recent paper 1 (hereafter referred to as I) we developed a scattering matrix analysis of Landauer conductance^ 
and of localization and reflection in coupled two- and three chain systems for weak disorder, using the Anderson 
tight-binding model. The full quasi-lD scattering set up consists, as usual, of the system of coupled disordered chains 
of length L connected to electron reservoirs by semi-infinite leads composed of non-disordered coupled chains attached 
at both ends. The leads carry the current which is incident on the disordered sample, together with the currents 
which are reflected or transmitted by the sample in modes of properly defined quantum channels 1 at the Fermi level. 

The analytical treatment of I was carried out for the case where the states at the Fermi energy in the various 
channels available for scattering of the injected electrons all belong to subbands of propagating states (i.e. the case 
where all channels are propagating). Thus, in this case all reflected- or transmitted electrons carry current. On 
the other hand, there may exist Fermi energy domains where the states at the Fermi level corresponding to incident 
electrons in a particular propagating subband include evanescent state solutions of the Schrodinger equations for some 
channels (evanescent channels) besides propagating solutions for others. 

In practice, Fermi energy domains such that propagating states in some channels coexist at the Fermi level with 
evanescent states in other channels are found when the interchain coupling exceeds the width of the subbands of 
propagating states in the leads. Evanescent states do not carry current, as is well-known*. They may however 
influence the current transport indirectly via their coupling, induced by the disorder, to propagating states at the 
Fermi level in other channels. 

The disorder-induced indirect effect of an evanescent channel in the Landauer conductance and in the related 
localization length 4 is quite distinct from the direct effect due to the disappearance of a conducting channel when the 
Fermi level is moved accross the edge of its propagating subband into an evanescent domain. Indeed, this process 
completely suppresses the primary (direct) effect which existed at Fermi energies within the propagating subband of 
the considered channel, where the latter was on the same footing as the other propagating channels coupled to it by 
the disorder. 

The object of this paper is to complement the discussion in I by a detailed analytical study of Landauer conductance 
and of localization lengths, as well as of reflection coefficients in Fermi energy domains where the relevant scattering 
modes in some channels are evanescent or non-propagating. As discussed above evanescent modes generally affect the 
conductance both directly, in that their existence implies the absence of corresponding (current carrying) propagating 
modes, and indirectly via their coupling to the other propagating modes. 

To our knowledge very few studies of the effect of evanescent modes in the conductance or the resistance of disordered 
wires have been published, for models which are different from the Anderson tight-binding model. Bagwell 5 studied 
in detail current transmission amplitudes and electrical conductance as a function of Fermi energy for electrons 
scattered from a single 5-function defect (and also numerically, for scattering from a finite range scatterer) in an 
otherwise non-disordered multichannel wire. On the other hand, Cahay et al& analysed the effect of evanescent states 
on the resistance of the two-dimensional random array of elastic scatterers numerically, using the scattering-matrix 
formalism of Datta, Cahay and Mc Lennani. 
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The scarcity of studies of the effects of evanescent modes (and the fact that the subject is not mentioned in any of 
the recent reviews or monographs discussing multichannel systems) does not distract from their intrinsic interest. It 
suffices to recall that the leads in the quantum conductance problem act as electron waveguides which define a basis 
for the scattering matrix of the multichannel disordered region. This indicates a similarity between the conductance 
problem and the study of optical waveguides where e.g. it is known that both propagating and evanescent wave 
solutions of the electomagnetic wave equation must be superposed in order to correctly describe the fields near 
sources or obstacles in an otherwise perfect waveguide 8 . On the other hand, in view of the generality of the Anderson 
model in the context of localization and transport in disordered systems, the present study of the role of evanescent 
modes should be of particular interest. Also, besides their relevance in various experimental situations^ few channel 
tight-binding systems are interesting because they can be discussed analytically for weak disorder, with the same 
degree of accuracy as corresponding ID-systems. 

In Sec. II. A we recall the Schrodinger tight-binding equations for the two- and three coupled chain systems studied 
in I. We briefly discuss the diagonalization of the interchain coupling terms carried out in I, which leads to the 
description of the leads in terms of independent channels for scattering of plane wave- and evanescent wave modes. In 
the three chain case we distinguish between equidistant chains on a planar strip with free boundary conditions and a 
system of equidistant coupled linear chains on a cylindrical surface. In Sec. II. B we recall the basic formulae for the 
Landauer two-point conductance and for the Lyapounov exponent of the conductance (inverse localization length^). 
These expressions are similar to those used in I except that the summations over channels in which an electron wave 
may be transmitted (reflected) are restricted to the only propagating channels. In Sec. Ill we summarize the main 
points of the determination of the transfer- and scattering matrices in these models referring to li for the detailed 
forms of these matrices (with proper adaptations for the case where propagating as well as evanescent channels are 
present at the Fermi level) . The final analytic expressions for the averaged coefficients of transmission- and reflection 
in propagating channels as well as expressions for localization lengths are included in Sec. IV. The discussion of these 
results together with further remarks is given in Sec. V. 



II. COUPLED TWO- AND THREE CHAIN DISORDERED WIRES 



A. Tight-binding models in channel bases 



Our A-chain model of a wire consists of parallel linear chains of Nl disordered sites each (of spacing a — 1 and 
length L = Nlcl) connected at both sides to semi-infinite non-disordered A-chain leads. 

The coupled two-chain Anderson model (A — 2) is defined by the Schrodinger equation in matrix form 



(1) 



where the ip l m denote the wave-function amplitudes at sites m on the ith chain, h is a constant matrix element for 
an electron to hop transversally between a site n on chain 1 and the nearest-neighbour site n on chain 2. The site 
energies £i m are random variables associated with the sites 1 < m < Nl of the disordered chain i, and £j TO = on 
the semi- infinite ideal chains defined by the sites m > Nl and m < 1, respectively. The above energies, including E, 
are measured in units of the constant hopping rate along the individual chains. 

The coupled three chain model (A = 3) is defined in a similar way by a set of tight-binding equations, whose 
actual form depends, however, on interchain boundary conditions. For free boundary conditions which correspond to 
arranging the parallel equidistant chains with nearest-neighbour interchain coupling on a planar strip, the Schrodinger 
equation is 



Vn+i + ^-A (E~e ln -h \ /VJA 

V&i + ti-i = -h E-e 2n -h \\<£\ , (2) 

with the sites in the disordered sections of length L = N^a and in the semi-infinite ideal chain sections labelled in the 
same way as in the two-chain case. On the other hand, in the case of periodic boundary condition which correspond 
to equidistant linear chains on a cylindrical surface, the tight-binding equations are 



<pl+i + ¥>i-i\ (E-exn -h -h \ / <pj\ 

^+1 + ^-1 = - h E-e 2n -h \\<&\ ■ (3) 
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An important aspect of the scattering matrix analysis of conductance in the above systems is the definition of 
bases for the scattering matrix which correspond to independent channels at the Fermi energy E. For this purpose 
one first transforms the Schrodinger equations (lll.'it to bases in which the interchain hopping terms are diagonalized. 
This transformation, discussed in I, replaces the coupled tight-binding equations describing the leads by decoupled 
equations corresponding to independent chain systems or channels. The transformed Schrodinger equations for the 
leads in the various systems under consideration arei 



Tpi+i + Tpl-A _(E-h o \ (ft 

ft+i + ft-i) - I E + h) [ft 



for the two-chain system 



(4) 



K ft +1 + 1 

for the three-chain system with the free boundary conditions, 



f ft n+x +ft n ^ 

ft+i+ft-i I = I E I I ft I , (5) 






n < 1 or n> N L , (6) 
for the three-chain system with periodic boundary conditions. The transformations U^ 1 to the new amplitude vectors 
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are given explicitely in I, as are the transformed Schrodinger equations for the disordered regions (equations (8.b-10.b) 
of I) in which the disordered channels are now coupled by the disorder. 

The bases for the definition of transfer matrices and the corresponding scattering matrices in the following section 
are provided by the Bloch wave (propagating)- and evanescent mode solutions of the Schrodinger eqations (|4I6|I of 
the form 



ft^ft^-e^ . (8) 

The propagating solutions for the various channels j correspond to real wavenumbers kj and exist in energy subbands 
defined by 



2 cos ki—E — h , 

2cosfc 2 = E + h , (9) 



for the two-channel systems, 



2 cos h = E - V2h , 
2 cos k 2 = E , 

2 cos k 3 = E + V2h , (10) 
for the three-channel system with free boundary conditions, and, finally, 



2coskx=E-2h , 

2 cos k 2 — 2 cos k% = E + h 



(11) 
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for the periodic three-channel system. On the other hand, the evanescent mode solutions for the various channels 
correspond to imaginary wavenumbers 

kj = iKj , cos kj = cosh/tj , (12) 
in ijHJ) and in (ftjll 1|> i.e. to exponentially growing or decaying solutions, 

^=<±~e=F*"* . (13) 

The evanescent modes in a given channel correspond to energies lying outside the energy band of propagating states 
in this channel. 

B. Conductance and localization 

We describe the conductance of a multichannel disordered wire by the Landauer two-probe formula^ 

g=^Tr(ii+) . (14) 

Here t denotes the transmission matrix associated with the M < N propagating channels of the iV-channel system, 
at the considered Fermi energy, 



(15) 

where Uj is the amplitude transmitted in a propagating channel i at the Fermi level at one end of the wire when 
there is an incident amplitude in the ( propagating) channel j at the other end. As in I the localization length is 
determined from the rate of exponential decay of the conductance^, namely 




J- = ~ 1™ Tr( ln 5) , (16) 

where averaging over disorder may be used, as usual, because of the self- averaging property of In g. 

The reason for restricting to the propagating channels in the double sum Tr(tt + ) in (|14|l is that evanescent modes 
at the Fermi level do not contribute to the current. However, since the conductance of the wire is defined as the 
current which flows through it divided by the electro-chemical potential difference between the reservoirs connected 
to the disordered region by ideal leads, the evanescent modes may play an indirect role through their effect on the 
electrochemical potentials&iS. Such effects have been discussed by Bagwell^ for electrons scattered from a single defect 
in an ideal quasi-lD wire. Here we consider similar effects for scattering from coupled chains of random atomic sites. 



III. SCATTERING MATRIX ANALYSIS OF TRANSPORT 



The analysis of transmission and reflection by the disordered region in I focused on the special case where all the 
scattering channels at the Fermi energy are assumed to be propagating. This situation is encountered typically when 
the interchain hopping parameter is less than half the width of the band of propagating states in the various channels 
i.e. \h\ < 2. In the opposite case, \h\ > 2, both evanescent and propagating channels will generally exist at the Fermi 
level. 

At the level of the general formalism only minor modifications of the analysis of I are required in order to incorporate 
the case where evanescent scattering channels may exist. For discussing this case we shall thus closely follow the 
treatment of I and refer to the latter for most of the details. The first step of the analysis involves defining transfer 
matrices for the disordered region of the quasi-lD systems described by (8.b-10.b) of I. These transfer matrices are 
used in a second step for obtaining scattering matrices giving transmission- and reflection amplitudes which allow us 
to study the conductance and related properties. 
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A. Transfer matrices 



Transfer matrices Y n for thin slices enclosing only a single site n per channel of the systems described by (8.b-10.b) 
of I are defined by rewriting these equations in the form 
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^-1 



v ; / 



(17) 



V ; ) 



The matrix Y n = Xq„ of dimension 4 for the case N = 2 and the matrices Y n = X' n , X n " of dimension 6 for the cases 
N = 3 with free boundary conditions, respectively, are given explicitely in equations (14-16) of I. 

The next step consists in transforming the Y n matrices in bases constituted in general by Bloch wave solutions 
(equation JSJ) of (1 IKill for some channels and evanescent wave solutions l|13|l for the remaining ones at the Fermi level. 
In these bases defined by transformation matrices W in I, the transfer matrices for the leads denoted, respectively by 
Yb = -^oo , Xq, Xq" for the considered cases are diagonal and of the forms 1 



Y = dmg(e lkl ,e- lk \e lk \e- ik \...) , (18) 

where the kj's, either real or purely imaginary (kj = iKj), are defined by l|9lll|) . Also, like in I, real wavenumbers kj 
are chosen to be positive, < kj < tt, so that e tkjn corresponds to a plane wave moving from left to right. Similarly a 
solutions ij} 3 n = e~ Kjn (Kj > 0) is viewed as an evanescent mode evolving from left to right (and, conversely, e Kjn as a 
mode evolving from right to left when discussing scattering from the disordered region. The full transfer matrices for 
thin disordered slices (acting on the column vector W -1 {ip„, V'n-i) ^ni ^n-ij ■ ■ •})) m the above mixed Bloch- and 
evanescent wave bases, denoted by Y n = Xo n ,X^ and X n " for N = 2 and for the two N = 3 cases, respectively, are 
given in equations (22) and (23) of I, in terms of tight-binding parameters defined by (22. a), (24) and (25) of I. These 
expressions remain indeed valid for Fermi energies corresponding to real values of some of the wavenumbers in l|9lll() 
and to pure imaginary values for others. 

Finally, the transfer matrices for the disordered wires of lengths L = N^a are products of transfer matrices in the 
mixed Bloch wave-evanescent wave bases associated with the Nl individual thin slices, 



Y L = [] Y n ■ (19) 
n=l 

As in I, the atomic site energies in (|1I3[I are assumed to be independent gaussian random variables with zero mean 
and correlation 



For weak disorder it is thus sufficient to explicitate (|19fl to linear order in the site-energies in order to study averages 
to lowest order in the correlation J5UJ. The latter implies indeed that different slices in l|19|) are uncorrelated. Under 
some notational proviso 11 - the final transfer matrices are given by (30) (N — 2) and (32) (N — 3) of I for the real or 
imaginary wavenumbers k\ , &2, k% in (|9I11|I at an arbitrarily chosen Fermi energy. 



B. Scattering matrices 



The scattering of plane waves (reflection and transmission) at and between the two ends of the random quasi- ID 
systems is governed by the 5-matrix, 



§ ~ (}++ 



(21) 
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where 

l 21 ^.2 ■ ■ ' , (22) 

and 



(23) 

Here tfj + (t~j~) and r ij + { r ij~) denote the transmitted and reflected amplitudes in channel i (which may be either 
propagating or evanescent) when there is a unit flux incident from the left (right) in a current carrying channel j. 
Left to right- and right to left directions are labelled + and -, respectively. The S- matrix expresses outgoing wave 
amplitudes in terms of ingoing ones on either side of the quasi-lD disordered wire via the scattering relations 



I'll 1 12 
'21 f o 



' 22 




Here I and /' (0 and 0') denote ingoing (outgoing) amplitudes at the left and right sides of the disordered region, 
respectively. It follows from current conservation that e.g. for a unit flux which is incident from the right in channel 
i one has£ 



£(l *JT I 2 + I rjt | 2 ) = 1 , (25) 

where the summation is restricted to the M < N current carrying channels at the Fermi level. Likewise, one also has 
for propagating channels j, 



M 

£(l t% + | 2 + | t%~ | 2 ) = 1 (26) 
i=i 

Note that in the case where channels are propagating the current conservation implies that the scattering ma- 
trix (|21|l is unitary [SS + = 1). This is no longer true, of course, in the presence of evanescent channels. As 
was shown in I the components of the out- and ingoing waves column vectors in (|24|l are given by the quanti- 
ties aj~ , a^ Q , • ■ • , a^ , a^ NL ,a+ Nh , a% NL and a^ , a% , a% , a^ NL ,a^ NL ,. . . , o,^ Nl , respectively, defined by 
components of wave transfer column vectors 

and obey the transfer matrix equation for the wires, 



{a liNL ,a liN ^,a 2iNL ,...} = Y L {a lfi ,a 10 ,a 2!0 ,...} . (28) 

Here the in- and outgoing waves include, in particular, evanescent waves evolving in directions towards- or 
away from the disordered wires. The components of the out- and ingoing waves column vectors in (|24|l are thus 
a i~oi a 2~o> ■ • ■ a N 0' a i L> a 2 Li ■ ■ ■ a ~N l an d a i~0' a 2 0i ■ • • a N 0' a i Li a 2 Li • • ■ a N Li respectively. With the so defined vec- 
tors of outgoing and incoming amplitudes, the S'-matrix is obtained by rearranging the equation (|28|) so as to bring 
them in the form l|24|l . The details of this somewhat lengthy calculation are explicitated in I. The explicit forms of 
the scattering matrices, Eqs (46-47) and (48, 48.a-48.f) of I, for N = 2 and N — 3, respectively, are expressed in 
terms of transfer matrix elements which are themselves defined in terms of general parameters given by (22. a), (24) 
and (25) of I with wavenumbers, real or pure imaginary, defined by (|9I11|) above. These S'-matrices readily yield the 
transmission and reflection submatrices in l|21|l . 
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IV. RESULTS 



In the appendix of I we obtained the explicit expressions of transmission- and reflection coefficients, \t- | 2 and 
|f*jj + | 2 , in terms of tight-binding parameters defining the transfer matrices of the disordered wires for cases where 
all channels are propagating at the Fermi level. In the appendix of the present paper we discuss the analogous 
expressions in more general cases where, some of the channels are conducting while the others are evanescent at the 
Fermi level. This requires explicitating the transmission- and reflection amplitude coefficients (|^ _ | 2 and |tyJ + | 2 ) for 
arbitrary wavenumbers, real or imaginary, for which the tight-binding parameters in (22. a) and (24-25) of I may thus 
be complex. In particular, the coefficients |^ - | 2 and |?"Vj + | 2 for conducting channels, which are relevant for studying 
the conductance and related transport properties, are influenced by the evanescent channels via the disorder mediated 
coupling between the two types of channels. 



A. Two-channel wires 



The two bands of propagating states in for real ki, and fc 2 are non-overlapping for \h\ > 2. In this case, to 
a propagating state in channel 1 corresponds an evanescent state at the same energy in channel 2 and vice versa. 
For the two-channel system the transmission- and reflection coefficients |^~| 2 and |f,^ + | 2 are given explicitely in 
1|A.1IA.3|) in the appendix, for the case where channel 1 is propagating and channel 2 is evanescent. By expanding 
these expressions to second order in the tight-binding parameters (|A.3|) and averaging over the disorder, using H20J1 . 
we obtain 

(fc-| 2 ) = i-^flr > ( 29 ) 

8 sin k\ 



8 sin fci 



(30) 



which obey the current conservation relation (|25|l for the conducting channel. 

From H29|) . (|14|l and (|16(l it then follows that the weak disorder localization length is 



L c 16 sin k\ 
The above results will be discussed further in Sec. V. 



1 h 

--—^TT- ■ (31) 



B. Three-channel wires 



Before discussing coefficients of reflection and transmission of plane waves in various propagating channels when 
evanescent modes are present, we first identify the various cases where both types of modes exist in the leads at the 
Fermi level, for free-boundary- and periodic systems, respectively. 

Free boundary conditions: if the subbands of propagating states in (|10|l for real ki,ks, k^ do note overlap, that is if 



\h\ > 2V2 , (32) 

than to an incident electron in a propagating state of energy E in a particular channel correspond evanescent states 
at the same energy in the other two channels. On the other hand, if the two outer bands 1 and 3 in 11U|) do not 
overlap while the inner band overlaps with the outer bands i.e. for 



V2<\h\< 2V2 , (33) 

there will be an evanescent mode present at the Fermi level for energies E lying within the overlap regions of the 
inner and outer bands. 
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Periodic boundary conditions: if the propagating band of channel 1 and the two degenerate bands of channels 2 
and 3 in lllfl do note overlap i.e. for 



, (34) 

then the waves at the Fermi level in channels 2 and 3 which correspond to an incident wave in the propagating band 
1 are evanescent. Conversely, for incident wave energies lying within the degenerate bands 2 and 3 the modes of the 
same energies in channel 1 will be evanescent. The same is true for \h\ < | for incident waves of Fermi energy within 
the non-overlapping parts of the band of channel 1 and the bands of channels 2, 3, respectively. 

In the following we will discuss successively the cases where one or two evanescent modes are present at the 
Fermi level for the two types of three-channel systems. For convenience the expression of transmission and reflection 
coefficients in terms of transfer matrix elements of the disordered regions given in (32) of I (under the proviso o£ii) 
is discussed in the appendix and diplayed in detail for the intrachannel transmission coefficients {\tjj~\ 2 jj = 1, 2, 3 in 
<IA.10IA.12J1. 



1. Free boundary condition 



a. Two evanescent modes 

Because of the symmetric arrangement of the propagating bands about E = we distinguish the cases where the 
Fermi level lies in the central fe-band or within one of the outer bands. 
For energies within the fe-band, we obtain from HA.11|) . 



and from (48.a), (48.e) and (32) of I, 



(&l 2 ) = i-^flr > ( 35 ) 

8 sin K2 



d^TH T^TT- - (36) 



i sin fc 2 



using l|2T?j) and (|A.8|I . It then follows from lfT^|l that 



1 



L c 16 sin &2 



The current conservation property l(2"K|l for propagating channels is clearly obeyed. 
On the other hand, for Fermi energies in the fci-band we get 



32 sin ki 

using (HJ and ifPJl . and from (48.a), (48.e) and (23) of I, 



-+|2\ 



M il I / ri • 2 ; ' 

32 sin fci 

which shows that (|25|l is obeyed. In this case we obtain for the localization length 



(37) 



(39) 



1 V 2 

4- = % ■ (40) 
L c 64sin 2 fci v ; 



b. One evanescent mode 
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As a typical case we choose the Fermi level within the overlap region of the fci- and fc2-bands, so that the matching 
mode in channel 3 is evanescent i.e. k 3 — in%. In this case we obtain successively from (|A.10() and (|A.11() . using H2UI) 
and l|A.8|) . 



32 sin ki \ sin ki sin fc 2 



! sin &2 V sm ^2 sin k\ 

and from (48.a-48.b) and (48.e-48.f) and (32) of I 



and 



(fe-| 2 ) = (I^T| 2 )= 1fi /f° , , (43) 
lb sm Ki sm k 2 



3Nl£ ° (44) 



32 sin k\ 

An important check on the correctness of equations (|41I46|) is to note that, again, they verify the fundamental sum 
rules i|25|) and (25. a) for the currents. The localization length for weak disorder obtained from the above results is 



1 



L c 32 \ 4 sin fci sin fc 2 smfcisinfo 



2. Periodic boundary conditions 



(47) 



For periodic boundary conditions, the current conservation properties (|25l -25.a') can be obeyed only if the disorder 
is restricted to identical realizations in the chains 1 and 2 i.e. 

ein = eani n = l,2...N L , (48) 

while the random site energies on chain 3 remain independent. This was shown in detail in I, in the case where all 
three channels arc conducting but remains true as well when evanescent channels are present, as expected. 

For energies in the propagating band of channel 1 in cases where there is no overlap with the degenerate propagating 
bands 2 and 3 (or for energies in the non-overlapping part of band 1 when this band partly overlaps the bands 2 and 
3) one has evanescent states at the Fermi level in channels 2 and 3. In this case we obtain, using H2U|) and (| A. 9|) . in 
(lA~10ll . 



rn 2 ) = i-^f?- , (49) 



and in (48.a), (48.e) and (32) of I, 



36 sin 2 k\ 



5NL f° . (50) 
36 sin fci 



10 



By inserting ij49jl in ij!6j) and expanding to linear order in the correlation Eq we find for the localization length 

1 5>=- 2 

4- = A ■ (51) 
L c 72sin 2 fci 

On the other hand, for Fermi energies such that the degenerate channels 2 and 3 are conducting and channel 1 is 
evanescent we find, from I|A.11IA.12() . 

<&l 2 > = <I^Tl 2 > = i-T^r > ( 52 ) 

4 sm fco 



and from (48.c), (48.e) and (32) of I, 



feTl 2 ) = (I* 3 TI 2 ) = t^4t . 

18 sm K2 



using (J201 and (|A.9|1 . Furthermore, from (48. a), (48. e) and (32) of I together with (HOI and HA.9|> we get, 

{\r^ + \ 2 ) = {\r^\ 2 ) = ^4r > ( 54 ) 
3o sm K2 



<k 2 /| 2 ) = (k 32 f | 2 )-#4 r ■ (W) 
18 sm fc 2 



The current conservation relation 1251 -25. a) are obeyed again and by using the above results in (j!6j) we find 

L c 72sin 2 fc 2 

V. DISCUSSION AND CONCLUDING REMARKS 

A simple interpretation of the results of Sect. IV for reflection- and transmission coefficients in two- and three 
channel systems in the presence of evanescent channels results from the comparison with the results of I for the case 
where all channels are propagating at the Fermi level. Thus we find that equations 129|) and l|30|) for the two-channel 
systems (N = 2) where the channel 2 is evanescent, follow exactly from equation (52) and (55) of I by suppressing the 
effect of the channel 2 (or, equivalently, assuming this channel to be absent). Similarly, we find that equations l|35|) 
and Ij36|) for the = 3 case with free boundary conditions follow from (62) and (68) of I by ignoring the effects of 
channels 1 and 3, which are now assumed to be evanescent. In the same way (|38ll and (|39|) follow from (61) and (67) 
of I by suppressing the effects of channels 2 and 3. Finally, in the case of a single evanescent channel 3 the equations 
(|41I46|) for the free boundary model follow from (61-63), (67-68) and (70) of I, respectively, by ignoring the effect of 
the channel 3. Exactly similar conclusions are obtained by comparing the results of Sec. IV for the N — 3 case with 
periodic boundary conditions with the corresponding results of Sec. IV of I. 

As we now discuss these simple properties of the reflection and transmission coefficients in Sec. IV are partly the 
consequence of destructive quantum interference, for sufficiently weak disorder of the terms describing the disorder- 
mediated coupling between propagating- and evanescent modes at the Fermi level in the corresponding amplitudes 
of transmission and reflection. This is illustrated e.g. in the equations l jA.10IA.12j l for the random intrachannel 
transmission coefficients |^~| 2 : the final double sum in each one of these expressions involves couplings between the 
propagating channel j and at least one evanescent channel assumed to be present at the Fermi Level. Now, since the 
terms m = n (the only once which survive in the averaging with I j20jl 1 in the sums over coupling terms of propagating 
and evanescent modes are pure imaginary (see parameters in ljA.8j) and (jA.9jl for kj real and fcj imaginary), they 
add up to zero in the intensity coefficients ljA.10IA.12j) . In other words, these coupling terms, while existing in the 
transmission amplitudes (as they arise from corresponding terms in ljA.7jl ) interfer destructively as a result of their 
special phases. 
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On the other hand, it is clear that the absence of coupling effects between propagating- and evanescent modes in the 
results of Sec. IV is specific to our weak disorder approximation. Indeed corresponding perturbative coupling terms 
at higher orders in the parameters l|A.8(l and IjA.Qjl in the transmission- and reflection amplitudes would have different 
phases, yielding non vanishing effects in the intensity coefficients. Similar higher order coupling effects (of fourth 
order in the random potential) between a single propagating mode and evanescent modes for a narrow strip-shaped 
pure wire with a random boundary ("surface") potential have recently been discussed by Makarov and TarasovJi. 

The localization lengths when evanescent modes are included, namely equation l|31f) in the two-chain wire, equations 
(|37H . (|40|l and l|47|l for the three chain wire with free boundary conditions and, finally, equations l|51|l and i|56|) for 
the three chain case with periodic boundary conditions differ qualitatively from the results of I (equations (58), (73) 
and (86) ini, respectively) for the usually considered case where evanescent states are absent at the Fermi level. The 
numerical coefficients in the above expressions indicate that localization lengths in the presence of evanescent states 
are generally enhanced by factors which vary typically between 1 and 2, with respect to corresponding values in I 
when evanescent states are absent. Strongly enhanced localization lengths in the presence of evanescent modes have 
also been observed in the numerical calculations of Cahay et al& for a different model. 

The results in Sec. IV for averaged reflection coefficients in the presence of evanescent modes may be used for 
deriving mean free paths for elastic scattering of an electron, using the formula^i^ 



1 1 M 

where M is the number of propagating channels in the quasi-lD conductor. By inserting the averaged reflection 
coefficients for the various quasi-lD systems discussed in Sec. IV, we find in all cases that 



4 = ^ , (58) 

where L c stands for the corresponding localization lengths for weak disorder given by (|31|) , (|37|) , l|4U|) , 147|) , l|51|l and 
(|52|l . respectively. These i c -values determine j- in the Born approximation for scattering by the disorder. We note 
that for systems with a single propagating channel (|58|l coincides formally with the relation between the localization 
length and the mean free path in a one-dimensional chain obtained by Thoulessi^. On the other hand, for systems 
with two propagating channels (|58|) is analogous to the general formula L c — M£ c , of Thoulessi^ for a quasi- ID 
system with M propagating channels. The equation l|58|l shows that, for weak disorder, the enhancement of the 
elastic mean free path due to evanescent states at the Fermi level is proportional to the corresponding enhancement 
of the localization length. 

The results of Sec. IV also allow us to express the conditions for the validity of the weak disorder scattering analysis 
of transport in quasi-lD systems in a precise physical form. Clearly our treatment is valid for intrachannel transmission 
coefficients close to unity (near-transparency) and sufficiently low reflection- and interchannel transmission coefficients. 
Roughly speaking this requires L = N^a < L c (a condition which can be made more precise for the various systems 
and Fermi energy domains by combining the results for the (|i^~| 2 ) in Sec. IV and in I with the corresponding 
expressions for L c ), which corresponds to the weak localization or (quasi) metallic regime. As is well-known, this 
regime allows one to find the correct expression for the localization length for weak disorder (eq << 1). This has been 
demonstrated, in particular, for the Anderson model for a one-dimensional chain where the localization length has 
been calculated analytically both in the weak- and in the strong localization regime (where all states are localized on 
the scale of L i.e. L > L c ) for weak disorder: the expressions obtained for both cases are identical and coincide with 
the familiar Thouless formula, Eq. (51) of I. 

In the present paper and in I, we have determined averaged conductances which allowed us to find localization 
lengths. It would be interesting to generalize these analyses to study conductance fluctuations and the ubiquitous 
universal conductance fluctuations in the considered coupled Anderson chain systems. In particular, since conductance 
fluctuations involve quartic terms in the tight-binding parameters (|A.8IA.9|) they would clearly be influenced by the 
disorder mediated coupling between evanescent- and propagating modes mentioned earlier in this section. On the 
other hand, it would be interesting to study localization in many-chain tight-binding models of quasi-lD wires. The 
exact analytical results obtained in this paper and in I for two- and three chain systems would be useful limiting cases 
for future treatments of localization in many-channel tight-binding wires. 
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APPENDIX: CALCULATION OF TRANSMISSION AND REFLECTION COEFICIENTS 

As in I we confine ourselves to discussing transmission- and reflection amplitude coefficients t~~ and r~- + associated 
with outgoing waves at the left of the disordered region. 

1. Two-channel wires 

In the case where k\ is real and k 2 — * K 2 is pure imaginary we are interested in finding 

I f— |2_ I --+ |2_ i^2l[ 2 (AU 

I 11 I _ m2 ' I r ll I - |j|2 ' { A ' L > 

where Xij denotes the matrix elements of the transfer matrix of the disorder region defined by (30) in I and 8 — 
X 22 X^-X 2 ^X^ 2 . The generalization of (A. 3) of I yields 

N L N L 



\ 6 \2 e -2, 2 N L = 1 + 2Im g ^ + g \a lm a ln + a 2m a* 2n + {b m b n e-^+^ {m -^ + c.c.) 

jn=l m,n— 1 

\X M \ 2 e~ 2 ^ N - =1 + J2 fl 2™«2„ (A.2) 

m,n 

\X 2 i\ 2 = oimiin cos2(™ - n)ki , 



where 



_ £ln + S 2n Sin + e 2 

O-ln — ——. — ; , a 2 



4 sin fei ii sinh k 2 

S 2 n — Sin 



(A.3) 



\\Ji sin k\ sinh K2 
2. Three-channel wires 

The transmission- and reflection amplitudes of interest defined by the matrices Si and S3 in (48. a) and (48. e) of I 
yield 

I^--I 2 = ^f . (A.4) 

where the are given by quantities 0k defined in I: an — 0i,ai2 — /?4,«i3 = 07, 

"2i = 02, ol 22 = 5 , a 23 = 0s, a 3 i = /3 3 , a 32 = fa, "33 = 

To second order in the site energies the expression (|A.4|) reduce to 

lNL \ 2 I Y 2 i | 2 , I r r + | 2 = |e^| 2 I F 23 | 2 , I TT+ | 2 = |e^| 2 | Y 25 | 2 , 
2 ^| 2 I Yu | 2 , I r 2 -+ | 2 = |e^| 2 I y 43 | 2 , I r 2 - 3 + | 2 = |e*^| 2 | F 45 | 2 , 
3 ^| 2 I ^1 I 2 , I r 32 + | 2 = l^^l 2 I ^63 I 2 , I r 3 - 3 + | 2 = |e^| 2 \ Y 65 | 2 . (A.5) 

Here the are matrix elements of the transfer matrix (32) in I of the disordered region and the /3j's and A are 
quadratic and cubic forms in these elements defined in (48. e) and (48. f) of I. The exponential coefficients in i|A.5|) 
differ from unity for imaginary wavenumbers. 

In the three-channel case we have a variety of different domains of Fermi energies where besides one propagating 
channel at least, there exist different evanescent channels, both for free boundary conditions and for periodic boundary 
conditions. Therefore we wish to explicitate the expressions l|A.4fl and l|A.5|) to second order in forms valid for arbitrary 







rli + 


| 2 = 








I 2 - 






-31+ 


I 2 - 
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wavenumbers k\ , fo, /C3, real or pure imaginary, which implies that the tight- binding parameters involved in the transfer 
matrix (32) of I will be generally complex too. 

On the other hand, for weak disorder, we only require the explicit forms of the quantities in ljA.4IA.5fl to second order 
in the site energies or, equivalently, to second order in the tight-binding parameters (24-25) of I. These explicit forms 
follow trivially for the reflection coefficients in lfA.5l) and for the interchanncl transmission coefficients in (fA.4jl since 
the elements of the reflection matrix f h given by (48. a) of I as well as the off-diagonal elements of the transmission 
matrix t given by (48. c) are proportional, to lowest order, to off-diagonal elements of the transfer matrix (32) in 
I, which are linear in the site energies. Also, as indicated earlier the form (32) of I of this matrix remains valid if, 
even in the considered case of complex wavenumbers and corresponding complex tight-binding parameters in If24l25jl . 
the symbols O* , with O = s,- , Uj , v^j , u>ij , are taken to mean replacement of the exponent coefficients i = 



'-I in 

(31) of I by — i. For these reasons we refrain from further explicitating the reflection coefficients and the interchannel 
transmission coefficients in HA.4IA.5f) . 



We now turn to the discussion of the explicit forms of the intrachanncl transmission coefficients, \t ■ 



,3 = 1,2,3, 



which play an important role and whose evaluation to second order requires more effort. Using the definition (48. e) 
and the explicit expressions of transfer matrix elements in (32) of I we obtain successively 



/3 5 = e -*(fci+k3)iVi 
/3g = _ e -*(*i+*>)Afr 



i(m— n)(fea— fea) 



m n m,n 



(A.6) 



Similarly we evaluate the cubic form A defined in (48. f) and (45) of I to second order in the tight-binding parameters 
in (32) of I. This yields 



1 - »y^(aim + b 2m + a 3m ) - ^2 



(m hi +n, foo +„, , + o „ •(m-n)(fe-fci) , . J(m-n)(k 3 -k2) , „ f p i(m-n)(fe 2 -fei) ] 



(A.7) 



The atomic tight-binding parameters given in (24-25) of I are: 



hn 
dn Qn 



£\n + 2£2n + £3r. 

8 sin k\ 

£\n + £3 



■1 c n — fn 



4 sin ki 

\/2(£ln " £3n) 

8%/sin k 2 sin ^3 



9n = Pn 



__ £\n + 2£2k + £3n 

8 sin k% 

V^jsin - £3n) 

8\/sin k\ sin k 2 

£\n — %£2n + £311 



8 V sin k\ sin k% 

for free boundary conditions, with the wavenumbers, real or imaginary, defined by IflUf) and 



(A. 



Q>\n 

d n 

Pn = 



£ln + £2n + £3n _ 2£2n + £3n 

■ ; 0>3n 



6 sin fci 

2ein + Esn 
6 sin ^2 

£3n — £2r, 



6 sin ki 



1 fn 



6 sin &2 

£\n — £3n 

6\/ sin k\ sin &2 ' 

2ei„ — E2n — £3n 

6\/sin fci sin fc 2 



9n 



£2n — £3n 
6\/sin ki sin &2 



~£\n + 2£2n — £3n £3n — £ln 



6-v/sin/c! sin & 2 



6 sin ki 



(A.9) 
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for periodic boundary conditions where the wavenumbers are given by (|ll|l . Finally we insert (|A.6|I and (|A.7(1 in i|A.4|) 
for the intrachannel transmission coefficients |i7j-~| an d evaluate the resulting expressions to second order in the site 
energies. This yields the relatively simple final expressions 



l*rr \ 2 \ e ~^ NL \ 2 



= 1 - 2 Im 2J «im + E [ax m a* n - 2 Re (ai m ai„)] 

- E [(.9™Pne l(m - n)(fe3 - fel) + c m f n e^ m - n ^- k ^ + c.c] , (A.10) 



|<--| 2 | e -^^ | 2 = 1-2 Im ^ b 2m + ^ [6 2m ^„ - 2 Re (6 2m 6 2 „)] 

X! [(rfm9ne l(m - n)(fe3 - fe2) + c m f n e l( - m ~ n ^- k ^ + c.c] , (A.ll) 



m m,n 



\^ 3 -\ 2 \e- lk3NL \ 2 



= 1 - 2 Im a i3m + ^ [a 3m a3 n - 2 Re (<23ma 3 n)] 

- E [(s™^" 1 -"^ 3 ^ + rf m 9ne l ( m - n )( fe3 -' £2 )) + c.c] , (A.12) 

m.n 

where summations over m and n run independently from m = 1 to m = Ni, and from n = 1 to n = N^. 
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